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One-Dimensional Electron Liquid in an Antiferromagnetic Environment: 
Spin Gap from Magnetic Correlations 
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We study a one-dimensional electron liquid coupled by a weak spin-exchange interaction to an 

. . . antiferromagnetic spin-S ladder with n legs. A perturbative renormalization group analysis in the 

0^ I semiclassical limit reveals the opening of a spin gap, driven by the local magnetic correlations on the 

0^ . ladder. The effect, which we argue is present for any gapful ladder or gapless ladder with nS 2> 1, is 

O^ ' enhanced by the repulsive interaction among the conduction electrons but is insensitive to the sign 

of the spin exchange interaction with the ladder. Possible implications for the striped phases of the 

P3 ■ cuprates are discussed. 

^ ; PACS numbers: 74.20.Mn, 74.20.-z, 75.20.Hr 

' , The coexistence of conducting electrons and localized spins remains one of the most challenging problems of con- 

^-H ■ densed matter physics, as evidenced by the enormous effort put into the study of, say, the Kondo lattice or doped 

,, antiferromagnets B. The recently discovered striped phases in La2-a:Sr2:Cu04 and various other high- Tc cuprates 

1—1 ■ add a new twist to this class of problems. "Stripes" is the name for spontaneously formed domain walls across 

^ [ which the two-dimensional antiferromagnetic order in these materials changes sign, and along which the doped holes 

5-H • are concentrated. The stripes are slowly fluctuating structures and may locally be modeled as metallic wires - in fact, 

C^ Luttinger liquids B - embedded in an antiferromagnetic environment. As suggested by Emery, Kivelson and Zachar 

Q , pair-tunneling of holes between the stripes and the environment may produce an electronic spin gap favoring either 

a charge density wave or superconducting correlations. Josephson coupling between stripes is expected to suppress 

the charge density wave, paving the way for superconductivity. Suggestions have also been made that a spin gap in 

the striped phase may be identified with the "normal-state" pseudogap observed in the underdoped cuprates B . 

In this Letter we also consider a one-dimensional electron liquid in an antiferromagnetic Mott insulating environ- 

Q ■ ment, and here focus on the role of the spin-exchange interaction between itinerant and localized electrons. This 

C^ ■ problem belongs to the more general class of Luttinger liquids in active environments l6l , a topic of importance not 

'~~' only to the striped phases, but also to e.g. nanotube and Kondo chain physics 1^. It is important to realize that 

I ' in the present case spin and momentum conservation severely restrict the possible relevant interactions between the 

^ , electron liquid and its environment. In particular, since the Fermi momentum of the Luttinger liquid (away from 

^^ ' half-filling) is incommensurate with that of any low-energy excitation of the Mott insulator we can neglect as irrele- 

^^ vant terms which transfer single holes to the insulator. Pair hopping is still allowed and is favored when the spins in 

^^ ■ the environment have a tendency to form singlets, as may be the case when there is a large pre-existing spin gap in 

the environment [H. In addition, however, a spin exchange interaction is always present, and is expected to become 

dominant for smaller gaps, correlating with a smaller density of local spin singlets. This is the case we consider here. 

^v' ■ Treating the localized spins semiclassically, we exploit a path integral formalism to construct a low-energy effective 

■^ , action with a companion set of perturbative renormalization group (RG) equations. Their solution reveals the opening 

C^ ■ of an electronic spin gap on the stripes, driven by the magnetic correlations in the environment. Rather strikingly, 

Ch , the effect is enhanced by the repulsive electron-electron interaction, but is insensitive to whether the coupling to the 

I ' environment is ferro- or antiferromagnetic. Although our approach allows for a fully controlled calculation only for 

' O ' large values of the localized spins or - as we shall see - for sufficiently wide antiferromagnetic domains between the 

^ : stripes, we shall argue that our results are robust in the limit of narrow spin-1/2 domains, at least in the case when 

the environment is noncritical. 

As lattice model we take a Hubbard chain (representing a stripe) coupled to the first leg of a neighboring spin 
ladder (representing the environment) by a spin-exchange interaction: 
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H = ^[-i(4+l,<TCr,o- + h.C.) -t- UUr^crrir^^cr + J rcI^^CT „^,Cr^^, ■ Sr,l 



r=l 



Jh{2_^ ^r,j ■ S^r+l.i + SrJ • SrJ + l)] , Sn+1,j — Sr,nug + 1 — 0, Jh > ■ (1) 

J = l 

Here Cr,a is a conduction electron operator at site r with spin index cr, Ur^a- = c]. a-Cr,^ is a number operator, and Srj 
is the operator for the localized spin at the site with coordinates r along legs and j along rungs. The model can be 
extended to include a coupling to antiphase ladders on either side of the stripe. As long as these ladders are correlated 
and the stripe is away from half-filling this will only change the magnitude of the couplings in the effective theory - 
to be derived below - but will not result in any qualitative changes [pi. 



The model in (|l]) can be taken as an effective model of a local stripe phase in the cuprates, valid on length- and 
time scales set by the fluctuation dynamics of the stripes (which is expected to be much slower than the dynamics 
of charge carriers along the stripe). For the purpose of exploring whether a spin gap opens up or not, we can count 
on the stripe as being metallic, as assumed in (|l|), since for weak disorder (induced e.g. by the dopant potentials) 
localization effects set in at length scales much larger than any relevant spin gap length scale. We should point out 
that since our model has the presence of stripes already built into it, the model cannot describe the instability that 
triggers the striped phases. For this, one must turn to other approaches, as in |10| . 

Given the Hamiltonian in (yj) , its partition function can be expressed as a Euclidean path integral by using coherent 
spin states in the semiclassical (large-S") limit of the localized spins, i.e. taking Sr.j 
of unit length. This gives 



Sftrj , where J7 is a vector 



V[n]V[c]V[c^]e 



-S[n,c,c^ 



(2) 



with action 
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S = dr iS 
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^Berry[^r.]{T)] + ^ cl,^drCr,a- + H{ 



,sn). 
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The first term in (^ is a sum over Berry phases, coming from the overlap of the coherent spin states: ^Berry['^r,j (t)] = 
L du nrj{u,T) ■ [du^r,j{u,T) X 9r fJ^.j (u, t)] , whcrc rir.j(u = l,T) = flrjir), flrj{u = 0,T) = Constant, with u a 
dummy variable. The Hamiltonian term H{c' , c, Sfl) in (^) acts at time slice r and is obtained from (|l]) by substituting 
electron and spin operators by corresponding Grassmann fields {c'jc) and classical vectors STJ, respectively. For the 
purpose of formulating a low-energy theory we linearize the electron spectrum close to the Fermi points ikp, assuming 
that U ^ ep ^ 2i(l — cosakp), and set Cr,a = \/a/^{e~'^^''°''^4'Lcr{af) + e^^''°-'^%l]B.cr{ar)), with a the lattice spacing, 
ne the electron density, kp = n^-K j^a^ and i^L/Ru left/right moving chiral fields. 

We expect that short-range antiferromagnetic correlations are present on the ladder also at the quantum level, 
implying that the partition function at low energies is dominated by paths with 






[{-ir+'^l-lljS^nr+tr.,IS)] 



(4) 



where Ur ■ trj = and \n\ = |ri| = 1. Here n is the local Neel-order parameter field, while f./ S represents small 
fiuctuations of the local magnetization [Q. For this to be a viable description of the ladder we require that the 
coupling to the conduction electrons is small, i.e. \Jk\ <C Jhi and also assume that the antiferromagnetic correlation 
length along the legs is much greater than the width of the ladder, allowing for n to be taken constant along the 
rungs fill. 

We first consider the case of free electrons [U — 0), away from half- filling (ne 7^ 1). Taking the continuum limit of 
(^ and neglecting terms of higher than quadratic order in i and 9^n, one obtains the action 



S = I dxdr 
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2TriS2_^{—iy{—n ■ {drn x dxTi)) (n x drfi) ■ "S^^j 



+ ^Hl^'dr + l^Vpdx)^ + —{Jl + Jr) ■ f-1 
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"0 = (V'Lj^fl)^ a Dirac fermion with velocity vp — 2atsmakp, and with 
^0 = cr^, ^1 = cr*', ■(/; = i/j^'^'y^. The Gaussian integral over £ in the partition function of (pi) can be carried out by means 
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J , with Ly 



JhS^j (6 - Sn - Sin,^^ ) + Jh5ij±i and Uj = -i{n x dru) + 6ji^J 



1 f 
with spin currents Jl —: ^i^L^f^a-fj.'ipL 

a"" 
of the substitution £[ 

where we define J± = J {J ■ n)n with J = J p + J p. We have here used the identity J ■ £ = {J — {J ■ n)n) ■ £ to 
preserve the constraint n{x) 1. £j{x) in the substitution £ ^ £ , an observation crucial to the subsequent analysis of 
the problem. 
This gives 



S — SNLa + Soirac + Si 



(6) 



where 

Snl^ ^ TT j dxdr i-{drnf + c{d^nf\ + 27riS'^(-l)^— / dxdr n ■ [drU x d^n) , (7) 

Soirac = ^ / dxdr^i^^dr + YvFdx)ip , (8) 

Si^Y I "^"^^^ [jKCii{n X drn) ■ Jj_ - ^C2Ji_ ■ J±] ■ (9) 

Here SnL(j is a nonlinear cr model describing the ladder, with coupling g^^ = S{JH'nieg'^^.i L~j )^'^ and velocity 

c = ciS{JHnieg/(^ij i^-^))^/^, and with the topological term SttzS"^ ■(— 1)-'^ J dxdrn ■ {drU x i9a;n) = i9Q, where 

9 = 27rS'X]?=i ("l)"* is the topological angle and QeZ the winding number of the mapping n : S^ ^ S^. Note that 
the topological term is absent for even-leg ladders and also effectively for odd-leg ladders with integer spin, while 
for odd-leg ladders with half-odd- integer spin it is present with 6 effectively equal to tt [|i3|. We shall return to the 
implications of this below. The Dirac action Soirac in (H) represents the electrons on the stripe, coupled to the ladder 
by Si in (|), with Ci = Y.^ Ki and C2 = L'^ ■ 

What is the effect of the interaction S*/? In particular, we wish to explore whether it may open up a spin gap for the 
electrons on the stripe. For this purpose we shall treat the interaction Si by means of a perturbative RG approach, 
using a mean-field formulation of the local Neel-order parameter field n. Specifically, we will derive an effective action 
for the spin sector which is valid over distances over which the spin ladder is ordered. Within the limits of validity 
of this action we then integrate out the short wavelength degrees of freedom to obtain its RG flow, allowing us to 
address the question above. 

Thus, given a patch in Euclidean space-time supporting local Neel order, we take the n-field to be in a fixed (but 
arbitrary) direction n. Introducing a local coordinate system (x,y,z) with z in the direction of n, and using the 
operator identity JlJl = ^Jl ■ Jl , we obtain from (P) — dropping the rapidly fluctuating first term of Si [p3 — an 

effective interaction Si valid up to length-scales of the size of the ordered region, 

^i = -^J dxdri^iJL ■ Jl + Jr ■ Jr.) + 2(JE J^^ + JlJ^)) , (10) 

with coupling gj = aJ^C'2/TT ~ aJ^/AnJii- Note that the spin anisotropy of the induced interaction in (10) is a 
direct consequence of the local Neel order of the n-field. Also note that the coupling gj is quadratic in Jk and hence 
the same for ferro- and antiferromagnetic spin exchange between the stripe and the environment. 

Bosonizing the Dirac action (H), i.e. splitting it into a charge boson and a (level fc = 1) Wess-Zumino-Witten 
Swzw,k=i model for the spin degrees of freedom, we absorb the quadratic terms of (|o|) into Swzw,k=i via a Sugawara 
construction, thus obtaining an effective action Sspin for the spin sector of the conduction electrons: 

Sspin — SwZW,k=l + Aq / d x{Jf^J^ -t- J^J^) , (11) 

where a;° = VsT, Vs — vf — 2gj, and with dimensionless coupling Aq — —gj/irvs- Since |Ao| ^ 1 we can use standard 
perturbative RG techniques to analyze Sspin, and at one-loop level we arrive at the scaling equations 

= 2^A^A'= ,k^j^t, (12) 



rflnL 

for the couplings A* of the operators Jj^Jii, with L a short-distance cut-off. Using (|l^ ) to solve for the RG flow, we 
obtain the trajectories A^ — (A^)^ — Xq with A = A^ = A^, and thus the scaling equation for A: dX/dlnL — 2ttX{X^ — 
Aq)^''^, which upon integration gives arctan(y%A/Ao)^^^) = 27r|Ao| InL/o. Hence, |A| grows under renormalization 
and at the length scale where |A| ~ 0(1) the perturbative treatment breaks down. This scale - where the perturbation 
is of the same order of magnitude as the fixed point action and renders the theory non-critical - defines the correlation 
length ^s of the electron spin sector. Using |A(^s)| ^ 0(1) 3> |Ao| in the scaling equation for A we thus obtain 
S,s ~ ae^/'*!'^''!, with an associated spin gap 

A«!^e-i/4|Ao|, (13) 

a 

The formation of a gap in this model is confirmed by the fact that ( pT| ) corresponds to a fermionic low-energy formu- 
lation of a spin-i XXZ chain (with a U(l)xZ2 symmetry) fla. The growing couphng constant scenario corresponds 



to an Ising anisotropy Jz > 1 of the XXZ chain, for which the latter is known to have a Neel ordered groundstate 
with a broken Z2 symmetry and a mass gap. 

The procedure leading up to ( [l3| ) requires that the environment exhibits Neel order over length scales exceeding 
£,s. Here we have to distinguish between spin ladders described by (JT]) with vanishing topological term (even- leg and 
odd-leg ladders with integer spin) and those where the topological term is present with 6 — tt (odd-leg ladders with 
half-odd-integer spin). The behavior of the nonlinear a model without topological term is well established [|5|; it 
has a finite mass gap and is ordered over distances given by the corresponding correlation length S.a- In contrast, the 
behavior when 9 = n is not rigorously known, although the consensus is that the topological term drives a crossover to 
the critical k=l WZW model at a length scale also set by ^a |p6[- However, in the weak coupling regime the topological 
term is effectively inactive [|l^, and as a consequence there is no distinction between gapless and gapful ladders on 
length-scales shorter than ^g.. It follows that the condition ^^ < ^g. validating our analysis is the same for gapless 
and gapful ladders. Evaluating g we find g^^ « 0.36S'n/eg, which in the weak-coupling regime with £_„ ^ age'^'^'^ pq ] 
implies the consistency condition 

O.SSnieg > ^ » 1 • (14) 

While (|lj) shows that our perturbative RG calculation is well-controlled only for large spins or wide ladders, it is 
important to emphasize that the interaction Si in (^ is well-defined for any values of S or nieg. As the symmetry 
of Si does not change when tuning the values of S or nieg, we expect that the result for the spin gap in (|l3) is 
analytic in these parameters with corrections that remain subleading as long as no topological effects intervene. 
On the other hand, when 9 = it, a. violation of (|lj) may change the physics, as suggested by bosonization and 
DMRG results for the Heisenberg-Kondo lattice model {nieg — 1,S — 1/2) llj,^: No gap is found for ferromagnetic 
coupling pO[ while for antiferromagnetic coupling the combined gap for itinerant and localized electrons scales as 
~ exp(— const. {'kJh/'2:+vf/Jk))- It might be appropriate to add a note concerning the prospect that non-perturbative 
effects at length scales larger than S,a could possibly carry over to the electron liquid. Although we cannot rigorously 
exclude it, it seems improbable considering the fact that the spin sector of the electron liquid develops a mass at a 
length scale which is shorter than and independent of ^cr and as such the mass is already well established at the scale 
where non-perturbative effects from the ladder may come into play. 

Let us now include the electron-electron interaction in (n^) [U ^ 0). At the level of the effective action for the 
electron spin sector this changes Sspin in ( |ll|) into 

Sspin = Swzw,k=i + d x\q{J^J^ + J^^Jfj) + ^^JlJr ) (15) 

with renormalized velocity Vg = vp — '^g.j — gu and couplings Aq = —{g.j + 9u)/''^Vs, Ag = —gu/TrVg, where gu — aU jir. 
Carrying out the RG analysis as above we obtain the spin gap 

A^!^exp(-"/'~"'-^y°/'^^)), (16) 



where S\ — -^/Aq — (Ag)^. Thus, as shown in Fig. 1, a repulsive electron-electron interaction ([/ > 0) produces a 
larger gap, while for [/ < the outcome depends on the precise ratio between g^j and g ,. An interpretation of the 
surprising scenario of a decrease of the gap for [/ < due to the environment is that the competition between the 
attractive electron-electron interaction, which enhances on-site singlet pairing, and the Ising anisotropy (discussed 
above), which enhances local Neel order, frustrates the system and hence reduces the gap. It should, however, be 
noted that the actual vanishing of the gap at gu / g.j = — 1 cannot be rigorously concluded from our model as the 
self-consistency condition ^^ < ^^ in this case requires Sni^g — > 00. 
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0^ ^ ^ . . ^^^' ^ -i^ ^ gu [vf/5] 

FIG. 1. The spin gap A as a function of gu\ the solid line is for gj = uf/5 and the dashed line for gj = 0. 



In summary, we have shown that a one-dimensional electron liquid weakly coupled by a spin-exchange interaction to 
a spin ladder with Snieg ^ 1 develops a spin gap. The gap exhibits a strong dependence on the sign and magnitude 
of the itinerant electron-electron interaction, but is insensitive to whether the coupling to the ladder is ferro- or 
antiferromagnetic. A symmetry argument implies that these results hold for any gapful ladder or gapless ladder with 
Snieg ^ 1. Applied to the striped phases seen in the cuprates this may suggest that the local antiferromagnetic 
correlations in the insulating domains may conspire with the electron correlations on the stripes to produce a sizeable 
spin gap. Details and extensions will be published elsewhere. 
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